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Abstract 

An analysis of diffusion in a supercooled liquid based solely in the density of diffusive directions and the value of energy barriers 
shows how the potential energy landscape (PEL) approach is capable of explaining the a and /3 relaxations and the fragility 
of a glassy system. We find that the /3 relaxation is directly related to the search for diffusive directions. Our analysis shows 
how in strong liquids diffusion is mainly energy activated, and how in fragile liquids the diffusion is governed by the density 
of diffusive directions. We describe the fragile-to-strong crossover as a change in the topography of the PEL sampled by the 
system at a certain crossover temperature Tx . 
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1 Introduction 

The study of slow dynamics in disorder systems and, 
in particular, the study of the glass transition in super- 
cooled liquids is a topic of considerable interest in con- 
densed matter physics. In general, liquids are divided 
in two different classes depending on the properties of 
their glass transitions, strong and fragile [1]. Strong liq- 
uids experience a gentle increase in the relaxation upon 
cooling, often according to the Arrhenius law, close to 
the glass transition temperature Tg. On the other hand, 
fragile liquids, experiment a sharp rise of the viscosity, 
increasing by several orders of magnitude in a very nar- 
row interval of temperature close to Tg. The glass tran- 
sition temperature is not related to any dynamical tran- 
sition and is experimentally defined as the one where 
the value of the viscosity is 10^'^ P. For strong liquids, 
nothing special happens close to Tg and the glass tran- 
sition shows a conventional behavior where no transi- 
tion may be defined. However, fragile systems seem to 
present a divergence behavior close to Tg indicating that 
some kind of new dynamical mechanism may be respon- 
sible for the onset of the glassy phase. Theoretical at- 
tempts conducted to study this possible new transition 
have been mainly focused in two different approaches: 
the mode coupling theory (MCT) [2] , and the potential 
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energy landscape (PEL) [3] . 

MCT studies structural arrest in supercooled liquids as 
a purely dynamic singularity happening as a result of 
a feedback between shear-stress, diffusion and viscosity. 
The idealized MCT predicts structural arrest to take 
place at a temperature Tmc > Tg [4,5]. To restore ergod- 
icity below TjviC additional hopping or activated mecha- 
nism have been introduced into the theory [6] , avoiding 
the kinetic singularity. MCT accurately describes im- 
portant aspects of relaxation dynamics in liquids above 
their melting temperatures, in particular, the behavior 
of the intermediate scattering function and the a and /? 
relaxations. The f3 relaxation plateau is related to the 
time expended by the particle to break the cage formed 
by neighboring particles. The accuracy of the MCT pre- 
dictions above Tg has been verified experimentally and 
using computer simulations [7,8,9]. 

The potential energy landscape (PEL) approach has 
been studied using the concept of an inherent struc- 
ture (IS) [10] (i.e., the configurations at the local mini- 
mum of the system's potential energy) . Several numer- 
ical and theoretical studies have provided evidence of a 
relation between the dynamics of the supercooled liq- 
uid and the PEL. It was found that correlation func- 
tions display stretching in time in the same temperature 
range in which the systems explores local minima of the 
PEL with deeper and deeper energy [11,12], a thermo- 
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dynamic description of the supercooled liquid was per- 
formed in terms of the IS configurational entropy [13], 
fragility was related to properties of the PEL [14] and 
the diffusion process was analyzed in terms of the visited 
inherent structures [15]. 

There was no connection between both theories (MCT 
and PEL). Basically there has not been a clear definition 

of from the potential energy landscape point of view 
and a precise landscape-based definition of hopping and 
activated dynamics has been lacking [16]. This situation 
has recently changed due to the instantaneous normal 
mode approach to the PEL [17]. This approach relates 
the diffusive processes to the number of accessible paths 
in the multidimensional energy landscape [18,19,20]. In 
particular, a key point is the temperature dependence of 
the fraction of negative eigenvalues of the Hessian cal- 
culated at the saddle points of the PEL [21,22]. It was 
found that this fraction approaches zc;ro at Tmc, but 
is appreciable at larger temperatures. Since the nega- 
tive eigenvalues correspond to diffusive directions on the 
PEL, the following scenario was proposed for the dy- 
namics of supercooled liquids close to the glass transi- 
tion: For T > TyiQ the system lies close to saddles in 
the configuration space and the relevant dynamic pro- 
cess is a diffusion among the saddle points along paths 
at almost constant potential energy, where there is no 
need to overcome any energetic barrier ( "border dynam- 
ics" ) . That implies that the factors impeding free diffu- 
sion of the particles are related to finding these paths 
(entropic factors), rather than due to energetic factors. 
However, below T^ic, there arc few paths available and 
diffusion is dominated by hopping processes allowing the 
system to evolve from minimum to minimum (minimum- 
to- minimum dynamics). This implies that a sharp slow- 
ing down on the dynamics should take place close to Tj^c 
if the energetic barriers to be crossed by the system are 
high enough [23]. 

2 Analytic approach to the PEL 

To check if this proposed mechanism is correct we are go- 
ing to study the problem analytically, relating the prop- 
erties of the; PEL the density of diffusive directions and 
value of the energy barriers — directly to the fragile and 
strong characteristics of supercooled liquids, and to the 
a and /3 relaxations predicted by MCT. 

We consider particles following a Brownian motion in 
three-dimensional space. In the absence of any other kind 
of interaction the mean square displacement is given by 

{r\t)) = QT[t-{l-e-% (1) 

Here we set the Boltzmann constant fc = 1, and take the 
particle masses m = 1. For very short times <C 1), the 
particles behave as free particles (r^(t)) = 3Tt^, and for 



longer periods of time [t 1) the particles behave as 
diffusive particles in a random walk {r'^{t)) = 6Tt. 

For particles in supercooled liquids such as Lennard- 
Jones systems, silica or water, the evolution is much 
more complicated due to non-trivial interactions among 
particles. These non-trivial interactions produce a very 
complicated energy landscape, making analytical results 
very difficult to obtain. However, a lot of information 
about these energy landscapes has been obtained using 
powerful numerical techniques such as Molecular Dy- 
namics or Monte Carlo simulations. In particular, it has 
been proposed that the density of diffusive directions 
with temperature, k{T), follows a power law [21] 

kiT) = AiT-Tucr- (2) 

The typical value AE of the energy barriers in hopping 
processes close to the mode coupling critical tempera- 
ture, has been determined for different models [23,24]. 
In all cases it has been found that a good approximation 
is given by 

AE « IOTmc- (3) 

For long enough time (we consider a time t to be long 
enough if t > 1), the diffusion of a particle is no longer 
free, making the total system to evolve over the multi- 
dimensional PEL. From a thermodynamic point of view 
the diffusion of the particle takes place if the direction 
chosen on the PEL turns to be a diffusive one, or if the 
barrier is low enough to be overcome by means of an ac- 
tivated process. Considering Eqs. (2) and (3), the prob- 
ability of diffusion of the particle at a temperature T 
from a time t > 1 to a time t -|- 1 is given by 

^^diff = A{T - TucT + [1 - - rMc)'^]e-^^/^. (4) 

Considering PdifTi it is possible to relate the total value 
of the mean square displacement {B?{t)) (at time t > 1) 
with the value of the random walk diffusion (r^(t)) by 

{R\t)) = E ( • Z J) {r''ii))P'^i{l - Pmf-'. (5) 
If t < 1 the diffusion is the one given by a random walk, 

{R\t)) = {T\t)). 

3 Results for (i?2(t)) 

Next, we study {B^{t)) from Eq. (5) for four different 
cases. Our standard case is going to be a binary mixture 
Lennard-Jones (BMLJ) system with density p = 1.2, 
as the one studied in ref.[21] for which Tmc = 0.435, 
7 = 0.94 and A w 0.05. We will study four cases. 
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Fig. 1. \n{R^{t)) vs. \n(t) obtained from Eq. (5) for temper- 
atures (from top to bottom) T = 5, 4, 3, 2, 1, 0.8, 0.6, 0.55, 
0.5, 0.475, 0.466, 0.454, 0.444, and 0.439. The values of the 
parameters are 7 = 0.94, T^q = 0.435, AE = WTy^Q and 
(a)fc(T) = 1, (b)^ « 0.5, (c).4 w 0.05 and {d)A = 0. 



4 Results for D{T) 

Once {R^{t)) is known it is possible to determine the 
dilFusion coefficient D{T) as a function of temperature, 
considering that a straight line, with unit slope, fitted to 
the long time behavior of data in Fig. 1, intersects a ver- 
tical line at ln(t) = at a height ln(6r'). To study what 
is the mechanism in the PEL leading to the behavior of 
a supercooled liquid as strong or fragile, we study three 
different cases. 

• Case (a): We study a system with no diffusive di- 
rections, where the only possibility for the system to 
diffuse is an activated process to overcome the ener- 
getic barrier AE ~ IOT^q. We consider again a value 
^ = ^ A; = 0. 

• Case (b): We consider the opposite case, where no ac- 
tivation process is available (since the energy barriers 
are extremely high) and the only diffusion mechanism 
is the search for diffusive directions. In order to study a 
system like that we take AE ~ IO^Tmc and A « 0.05. 

• Case (c): Finally we consider the Binary Lennard 
Jones system considered in [21] where both diffusive 
mechanisms are available. 



• Case (a): We consider the simplest system where only 
Brownian motion is present and there is no effect of the 
PEL. We take k{T) = 1, implying that any direction 
chosen by the system is going to be a diffusive one. 

• Case (b): The second system considered has a larger 
density of diffusive directions that the one studied in 
[21]. To do so we take A « 0.5. In this particular case 
k{T) = 1 for T > 2.5. 

• Case (c): We consider the standard case from ref.[21] 
where k{T) ^ 1 for every T in this study. 

• Case (d): We consider a case where there are no dif- 
fusive directions A = k{T) = 0. 



Results are presented in Fig. 1 For temperatures rang- 
ing from T = 5 to T = 0.4386. A clear relation between 
the (3 relaxation time and the density of diffusive direc- 
tions in the system, k{T), is found in Fig. 1. In Case (a) 
all directions are diffusive and particles do not expend 
any /3 relaxation time searching for a diffusive direction 
to scape, that is the reason why there is no plateau in 
Fig. 1(a). On the contrary, in Case (d), there are no dif- 
fusive directions and the particle expends a long time 
in the (3 relaxation plateau. The only possible diffusion 
mechanism to scape from this plateau is by an activated 
process to overcome the AE' barrier. So the f3 relax- 
ation time may be interpreted in terms of the PEL as a 
"search" for diffusive directions. This mechanism should 
be equivalent to the one described in the MCT where 
the particles "search" for directions to scape from the 
cage formed by surrounding particles. 



Results are presented in Fig. 2. We plot ln(6D) vs. 1/T 
in Fig. 2(a) and ln(6D) vs. ln(T - Tmc) in Fig. 2(b). 
Note how Case (a) clearly shows a linear (Arrhenius be- 
havior) , typical of a strong liquid. Nothing special hap- 
pens when T = T^ic , since the mechanisms of diffusion 
are always the same (activated processes). On the con- 
trary Case (b) presents the typical behavior of a fragile 
liquid predicted by MCT. In this case, the transition to 
a glassy state at T = TyiQ, is a singular one. Since the 
only mechanism of diffusion are the diffusive directions 
and those are equal to zero at T = TyiQ the system 
gets trapped in a glass state by dynamical arrest. Note 
how the barrier must be very high to get a really slow- 
ing down in the dynamics, making any kind of hopping 
impossible. This result agrees with predictions made in 
[23] . If the energy barriers are not high enough we are in 
Case (c) and we have a transition from strong to frag- 
ile behavior at T = T^ic as the one predicted by [25] 
and observed numerically in [26]. Note how Case (c) is 
almost identical to Case (b) when the temperatures are 
far from Tj^q clearly indicating that the mechanisms 
of diffusion are governed mostly by "border dynamics," 
that is, by evolution of the system through the diffu- 
sive direction of the PEL. However for T < TyiQ, Case 
(c) behaves like Case (a), which means that the mech- 
anisms of diffusion are now governed by activated pro- 
cesses to overcome the AE energy barriers. Case (c) is 
a clear example of crossover from "border dynamics" to 
"minimum-to-minimum dynamics." 

Results in Fig. 2 are easy to understand considering 
that, since Pdiff is the probability for a particle to dif- 
fuse, a Pdiff 7^ turns t into a lower value, Pdiff*- That 
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Fig. 2. Figure 2(a) shows lri(6-D) vs. l/T. Dashed hnc is a 
hncr fit to the data for Case (a), dotted hue marks the value 
T = T^Q. Figure 2(b) shows ln(6i:') vs. ln{T-T^Q). Dashed 
line is a liner fit to the data for Case (b). The values of the 
parameters in the different cases are Case (a) AE = IOT^q, 
A = 0, Case (b) AE = IO^Tmc, A « 0.05 and Case (c) 
AE = WTyiQ, ^«0.05. 

means that instead of (i?^(i)) = 6Tt for t >> 1, we 
have (i?^(t)) = QTP^int. For the Case (a) (strong glass), 
since A = 0,we have Pdiff = e~^^/^, making (-R^(i)) = 
&Te~^^/^t, and getting a value of the diffusion coeffi- 
cient \Yi{D) = ln(r) - AS/T. For the Case (b) (frag- 
ile glass) it is possible to consider in good approxima- 
tion that Pdiif = A{T - TyicV, obtaining ln(D) = 
ln(T) + In(^) + 7ln(T - Tmc), which is the behavior 
expected by MCT. The constant value for AE given in 
Eq. (3) is valid only near Ty^Q. If we want to obtain re- 
sults from our model in a larger range of temperatures, 
a non-constant value of AE, as the one reported in [23] , 
should be taken into account. 



Equation (7) implies an Arrhcnius behavior for D(T), 
making impossible to find a crossover from fragile to 
strong. It means that there must be some crossover tem- 
perature Tx where k{T) changes from a power-law to 
an Arrhenius-Iaw, marking a change in the PEL topog- 
raphy and the beginning of a crossover from fragile to 
strong. Experimental results [30] have shown that Tx 
normally has a value very close to TyiQ. 

To study the effect of this possible change on topography 
we arc going to modify the density of diffusive directions 
of Case (C) , considering a new Case (D) with 

k{T) = /oe^^[(i/^MC)-(i/T)] T < Tx , 

k{T)= A{T-TMcr T>T^, (8) 

where /o is the density of diffusive directions at TyiQ, 
which is now set to a value different from zero (/o = 
0.001). The crossover temperature, Tx, is given by the 
lower root of the equation 

/oe^^[(V^MC)-(VTx)] = ^(r^ _ Ty^^y. (9) 

A plot of fc(T) vs. T is shown in Fig. 3a compared with 
the one corresponding to Case (C), where A;(Tmc) is 

strictly equal to zero. 

Results for ln(6T') vs. l/T are given in Fig. 3b and com- 
pared to the behavior of the pure fragile system. Now the 
crossover from fragile to strong takes place at T = Tx 
(close to Ty[Q) and it is not so abrupt as the one in 
Fig. 2a. The qualitative behavior shown in Fig. 3b coin- 
cides with the result found for Silica by means of Molec- 
ular Dynamics Simulations [26]. However, the behavior 
reported in Fig. 2a is more similar to the one found for 
BMLJ [24], posibly indicating that for BMLJ /o « and 
Tmc ~ ■ 

6 Conclusions 



5 Crossover from Fragile to Strong 

Numerical simulations and theoretical calculations have 
recently shown that the density of diffusive directions 
is not exactly zero at T^/[c [27,28]. It has been argued 
that this density behaves as an Arrhcnius exponential 
decay which is only zero at T = [29] . Actually, a close 
inspection of the data reported in Rcf. [26] shows clearly 
that k{T) is almost null at TmC) but not exactly zero. 

If we change the power law behavior in Eq. (2) by an 

exponential decay given by 

k{T) = A*e-^^'l^, (6) 
with A* constant and AE* an energy scale, we obtain 



To conclude, we have analyzed the diffusion in a super- 
cooled liquid based solely on the density of diffusive di- 
rections and activated processes and have shown how 
the PEL provides an explanation for the (3 and a relax- 
ations and the fragility characteristics of a glassy sys- 
tem. The [3 relaxation is directly related to the attempts 
of the system to move through the diffusive directions. 
The PEL shows that a "strong" liquid is one in which 
the main mechanisms of diffusion are "activated dynam- 
ics" and that a "fragile" liquid exhibits dynamics typical 
of sup{Tcx)oled liquids with diffusive directions but very 
high barriers where hopping is almost impossible. In this 
case PEL supports the same dynamical arrest behavior 
predicted by MCT. The crossover from fragile to strong 
is found to be related to a change on the topography of 
the PEL at a certain crossover temperature Tx , where 
the density of diffusive directions changes from power- 
law to Arrhenius. 
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Fig. 3. Figure 3(a) Density of diflFusive directions k(T) vs. 
temperature T for Case (C) (thin line) and Case (D) (thick 
hne). T^Q marks the value k{T) = for Case (c) and Tx 
marks the change on topography of the PEL for Case (D). 
Figure 3(b) shows ln(6£)) vs. 1/T for Case (D) (black points). 
Continuous line represents the behavior for the pure fragile 
glass former. Tj^q and Tx are marked with dotted lines. 
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